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We argue that the Square Kilometre Array has the potential to make both redshift (HI) surveys and radio 
continuum surveys that will revolutionize cosmological studies, provided that it has sufficient instantaneous field- 
of-view that these surveys can cover a hemisphere (/sky ~ 0.5) in a timescale ~ 1 yr. Adopting this assumption, 
we focus on two key experiments which will yield fundamental new measurements in cosmology, characterizing 
the properties of the mysterious dark energy which dominates the dynamics of today's Universe. Experiment I 
will map out ~ 10®(/aky/0.5) HI galaxies to redshift z ~ 1.5, providing the premier measurement of the clustering 
power spectrum of galaxies: accurately delineating the acoustic oscillations and the 'turnover'. Experiment II will 
quantify the cosmic shear distortion of ~ 10^" (/sky /0. 5) radio continuum sources, determining a precise power 
spectrum of the dark matter, and its growth as a function of cosmic epoch. We contrast the performance of 
the SKA in precision cosmology with that of other facilities which will, probably or possibly, be available on a 
similar timescale. We conclude that data from the SKA will yield transformational science as the direct result 
of four key features: (i) the immense cosmic volumes probed, exceeding future optical redshift surveys by more 
than an order of magnitude; (ii) well-controlled systematic effects such as the narrow 'fc-space window function' 
for Experiment I and the accurately-known 'point-spread function' (synthesized beam) for Experiment II; (iii) 
the ability to measure with high precision large-scale modes in the clustering power spectra, for which nuisance 
effects such as non-linear structure growth, peculiar velocities and 'galaxy bias' are minimised; and (iv) different 
degeneracies between key parameters to those which are inherent in the Cosmic Microwave Background. 



1. Background 

Over the last few years we have entered an 
'era of precision cosmology' in which rough esti- 
mates of many of the key cosmological parameters 
have been replaced by what can reasonably be de- 
scribed as 'measurements' (parameters known to 
better than ~ 10 per cent accuracy). Remarkable 
progress in observations of the temperature fluc- 
tuations in the Cosmic Microwave Background 
(CMB) has been central to this scientific transfor- 
mation (e.g. the recent results from the WMAP 
satellite, Spergel et al. 2003). 

The current standard cosmological model 
('ACDM') assumes that the Universe is composed 
of baryons, cold dark matter (CDM) and Ein- 
stein's cosmological constant A. The present-day 
densities of baryons, (baryons -I- CDM) and A, as 
a fraction of the critical density, are denoted by 
fib, ^^m and r^A, respectively. The local rate of 



cosmic expansion is described by the Hubble pa- 
rameter h = i7o/(100 km s^^ Mpc~^). Structures 
in the Universe are assumed to grow by gravita- 
tional amplification of a primordial 'power spec- 
trum of fluctuations' created by inflation, which 
is taken to be a 'scale free' power law Pprim(fc) oc 

"scalar 1 

Within this framework, we begin by reviewing 
what the CMB tells us about the cosmological 
model. The basic CMB observable is the angu- 
lar power spectrum of temperature fluctuations, 
which consists of a series of 'acoustic peaks'. We 
note that the physics of the CMB fluctuations (i.e. 



power law function is 'scale free' because it contains no 
preferred length scale. The power law index ng^alar = 1 
is described as 'scale invariant' because in this case, the 
matter distribution has the same degree of inhomogeneity 
on every resolution scale (is 'self-similar'). Inflationary 
theories of the early Universe predict that this should be 
approximately true. 
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the photon-baryon fluid driven by the dark mat- 
ter potential wells) depends entirely on the phys- 
ical densities of matter and of baryons, il^h"^ and 
Qbh^.^ However, the projection of the anisotropy 
spectrum onto the sky (i.e. the angular locations 
of the acoustic peaks) also involves the angular 
diameter distance to the last-scattering surface. 

• The angular location Oa of the first 'acous- 
tic peak' of the temperature power spec- 
trum turns out to be relatively insensitive 
to all parameters apart from the curva- 
ture of the Universe. The global geome- 
try is close to spatially flat, although pos- 
itive spatial curvatures « —0.05) and 
hence closed-Universe models remain con- 
sistent with all available data (Efstathiou 
2003). It is important to note that exact 
spatial flatness, typically assumed, remains 
a theoretical prior which demands observa- 
tional confirmation. 

• The physical matter density {Cl^h^ « 0.14) 

is set by the ratio of the first acoustic peak 
height to the large-scale amplitude of the 
temperature power spectrum. 

• The physical baryon density (fib^^ ~ 0.024) 
is fixed by the relative heights of the first 
and second acoustic peaks (and is nicely 
in accord with the value derived from 
nucleosynthesis-based arguments) . 

• The values of flmh"^ and fl^h^ together de- 
termine the sound horizon s at recombina- 
tion, the characteristic physical scale fix- 
ing the spatial position of the first acous- 
tic peak. Using the angular location 9^ 
of this peak, we can derive the angular di- 
ameter distance Da = s/6a to the surface 
of last scattering (whose redshift z 1100 
is known from thermal physics). Assum- 
ing a flat Universe and the ACDM model, 
Da{z = 1100) depends on a combination 

^The factor h'^ arises because these physical densities are 
proportional to the critical density of the Universe, which 

scales as H^. The existence today of non-zero curvature 
or of a cosmological constant A is insignificant at recom- 
bination. 



of and h. Together with the value of 
flmh^ known from the first peak height, this 
breaks the degeneracy between and h, 
yielding ~ 0.3 and h « 0.7. 

• The overall normalization of the tempera- 
ture power spectrum (i.e. the fluctuation 
strength at z = 1100) can be related to 
the amplitude of mass fluctuations at red- 
shift zero, conventionally defined as the rms 
mass perturbation in spheres of radius Sh^^ 
Mpc assuming linear evolution'' (and de- 
noted (Js). In a ACDM model, if the matter 
densities and Hubble constant are known, 
then this quantity is still degenerate with 
the optical depth to the last- scattering sur- 
face (r). The detection of polarization in 
the CMB by the WMAP experiment begins 
to break this degeneracy by obtaining the 
value r ~ 0.17 (hence as ~ 0.9). 

• Changing the slope of the primordial power 
spectrum of mass fiuctuations (ngcaiar) in- 
duces a broad-band tilt in the CMB temper- 
ature power spectrum, that may be quanti- 
fied by e.g. the ratio of the first and third 
acoustic peaks (note that the quantity flbh'^ 
only controls the ratio of odd-to-even peak 
heights). Current data favours a roughly 
scale-invariant primordial fluctuation spec- 
trum (riscalar « 1)- 

In summary, according to this crude argu- 
ment, the CMB yields six independent quanti- 
ties (assuming polarization measurements) . For a 
flat ACDM model, the WMAP experiment alone 
measured values of (rim, ^b, h, a^, r, rigcaiar) to 
5—10 per cent accuracy (see Spergel et al. 2003). 

Similar values for some of these parameters are 
yielded by several independent, non-CMB tech- 
niques, for example: 

^When the amplitudes of density fluctuations are small 
(when the overdensity 5 <SC 1), the equations of struc- 
ture formation may be linearized using perturbation the- 
ory. Gravitational growth is said to be in the linear regime 
and may be solved analytically with relative ease. Eventu- 
ally, the fluctuation amplitudes will become too high and 
linear theory will breaJj down. 
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• Observations of Cepheid variables and sec- 
ondary distance indicators to measure h 
(Freedman et al. 2001). 

• Searches for high-redshift supernovae to 
measure Ha (Riess et al. 1998; Perlmutter 
et al. 1999). 

• Measurements of the abundance of rich 
clusters to determine as (e.g. Bahcall & 
Bode 2003). 

• Measurement of the amplitude of cosmic 
shear (weak lensing due to large-scale struc- 
ture) to constrain as (e.g. Brown et al. 2003; 
Hoekstra, Yee & Gladders 2002; Jarvis et al. 
2003; Bacon et al. 2003; Pen et al. 2003). 

The consistency of these results with those ob- 
tained from the CMB demonstrates that - if 
the standard cosmological picture is correct - 
these methods are not obviously compromised 
by systematic errors at the current level of accu- 
racy. The next-generation CMB experiment, the 
Planck satellite, is planned to launch in 2007 and 
will further refine CMB measurements of these 
six parameters and second-order effects via (i) 
drastically-improved sensitivity on small scales 
and (ii) an enhanced polarization capability."* 

However, if we wish to relax any of the assump- 
tions of our most basic (flat ACDM) model then, 
as we will describe, the degeneracies inherent in 
the CMB become insuperable and extra informa- 
tion is required (Efstathiou & Bond 1999; Bri- 
dle et al. 2003). Just such a course of action is 
demanded by today's most pressing cosmological 
questions: 

• Is the current rate of cosmic expansion ac- 
celerating, as suggested by recent supernova 
measurements? 

• If so, does Einstein's cosmological constant 
(vacuum energy) exist, despite our inabil- 
ity to understand its value theoretically? 
Or, is the accelerating expansion driven by 

''See http: //www. astro . esa. estec .nl/~Plaiick| The 

temperature power spectrum measured by WMAP is al- 
ready 'cosmic variance limited' (i.e. cannot be improved) 
on large scales up to and including the first acoustic peak. 



a different form of dark energy, such as 
'quintessence' ? What are the properties of 
this dark energy? 

• Can competing models of inflation be dis- 
criminated by accurate measurements of 
the shape of the primordial power spectrum 
(e.g. by a 'running spectral index' instead of 

'^scalar )'^ 

• Can the tensor (gravitational wave) compo- 
nent be isolated? 

• Is the Universe exactly flat or does it possess 
a small spatial curvature Slj,? 

In this paper we focus our discussion upon the 
quest to characterize the mysterious dark energy, 
which drives the accelerating cosmic expansion 
and contributes the majority of the present-day 
energy density of the Universe. As an explana- 
tion of the dark energy, Einstein's cosmological 
constant A is consistent with current observa- 
tions but poses tremendous theoretical challenges. 
The 'natural' (quantum-mechanical) value of A 
is a factor 10^^° higher than that measured 
by astronomers. This severe difhculty has mo- 
tivated the development of alternative models, 
many featuring a dynamic component of dark en- 
ergy whose properties evolve with time, such as 
quintessence. These competing models are essen- 
tially untested and require experiments able to 
measure the properties of dark energy as a func- 
tion of cosmic epoch. The dark energy model is 
commonly characterized by its equation- of- state 
w{z) = P/ p relating its pressure P to its energy 
density p (in units where the speed-of-light c = 1). 
For Einstein's cosmological constant, w{z) = —1. 

The resolution of this problem is the current 
cosmological frontier, and has been widely identi- 
fied as having profound importance for physics as 
a whole. Needless to say, if there were a paradigm 
shift in the next decade that alleviated the dark 
energy problem, then the SKA experiments de- 
scribed below would be even more important in 
confirming and better understanding the new cos- 
mological model. 

Can the CMB inform us about dark energy? 
Although in the standard (ACDM) model dark 
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energy has a negligible effect at recombination, 
the CMB is in fact sensitive to the equation-of- 
state of dark energy because the angular position- 
ing of the acoustic peaks depends on the angular 
diameter distance to z = 1100, which is affected 
by w{z) (Huey et al. 1999). However, the efficacy 
of this probe of w{z) is limited because: 

• Using CMB data alone, there is an imme- 
diate degeneracy between w(z) and Sim to 
produce the same value of Da{z = 1100). 

• Even given accurate external information 
about rim and h, the quantity Da{z = 
1100) can still only constrain a weighted av- 
erage oiw{z), i.e. one degree of freedom. 

• The precision with which this one degree of 
freedom can be measured is further limited 
because w{z) only influences dynamics at 
low redshift, thus its influence upon Da{z — 
1100) is smah. 

• This lack of sensitivity implies that we can 
make small adjustments to the value of the 
sound horizon s at fixed by changing 
the mix of baryons and CDM (i.e. fib) to 
leave unchanged the angular locations of 
the CMB acoustic peaks, which depend on 
s/Da[z = 1100). 

This difficulty in measuring a non-Einstein 
dark energy equation-of-state 'w{z) using the 
CMB alone, even if assumed to be a constant 
w'cons 7^ —1, is illustrated by Figure^] in which we 
simulate the degeneracy between i7m and Wcons 
after Planck data has been analyzed. Even with 
the addition of a very accurate independent mea- 
surement of ft,, there remains a ±0.1 scatter in 
the value of ifcons- We note that changing Wcons 
from —1 to —1.1 alters the value of Da{z = 1100) 
by less than 1%. This is because the co- moving 
distance to a given redshift is determined by the 
integral 

^ r (c/gp) dz' 

Jo V^n^il + Z')3 + (1 - a„)(l + z')^('+'"™"=^ 



to which dark energy only contributes at very low 
redshifts (z<l). For fixed il^ and h, the sound 
horizon can be changed by a similar fraction by 
varying fli^h^ within its uncertainty range after 
the Planck experiment. 

The prospects for using the CMB to map out 
any variations in the dark energy equation-of- 
state with cosmic time, the key observation re- 
quired to discriminate quintessence cosmologies, 
are even worse. Even if the exact value of 
were known, the CMB primary anisotropics only 
measure a weighted average of w{z) (Saini, Pad- 
manabhan & Bridle 2003). For example, in a 
simple model 'w(z) — wq + wiz, we can al- 
ways scale wi with wq in such a way to keep 
Da{z — 1100) fixed and leave the power spec- 
trum of the primary anisotropics invariant (al- 
though there will be second-order changes, such 
as Sunyaev-Zeldovich decrements due to galaxy 
clusters, and the late-time Integrated Sachs- Wolfe 
effect). 

We conclude that, due to degeneracies intrinsic 
to the CMB, the outstanding questions of the era 
of precision cosmology demand additional experi- 
ments of comparable power, with different inher- 
ent parameter degeneracies. As discussed above, 
the CMB can only weakly constrain the dark en- 
ergy model. An accurate measurement of w{z) 
requires extremely precise cosmological data at 
lower redshifts, where dark energy significantly 
affects the dynamics of the Universe; for exam- 
ple, per-cent level determinations of the distances 
to redshift slices at z ~ l.'^ What cosmological 
probes can provide such precision? 

• Observations of distant Type la supernovae 
as 'standard candles' track the luminosity 
distance as a function of redshift and pro- 
vide an 'orthogonal constraint' to the CMB. 
Supernova surveys will be a powerful probe 
of the dark energy model, but we argue that 
it is not yet proven that the experimental 
systematics can be controlled to sufficient 
precision to disentangle the effects of dark 

^Beyond z fs; 1.5, dark energy is expected to have a negli- 
gible effect on dynamics. This assertion should be tested; 
but if true, distance determinations to z 3 have lim- 
ited power for constraining ui(z), and measurements of 
the Hubble constant at z « 3 have no such power. 
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Figure 1. Simulated joint constraints (68%, 95%) on the matter density f2m and a constant dark energy equation- 
of-state w{z) — Wcons, resulting from Planck observations of the CMB. We use TT, TE and EE simulated Planck 
data and only retain points with f > 50 to account for uncertainty in any dark energy perturbations, which 
may manifest themselves on very large scales and are strongly dependent on the exact dark energy model. The 
strong degeneracy in the (fimjWcons) plane arises because the CMB power spectra are approximately invariant 
for fixed values of (1) Qmh^ and (2) the angular diameter distance to the last-scattering surface, Da{z = 1100), 
which fixes the angular location of the first acoustic peak. Da{z — 1100) may be held constant by varying Qm 
simultaneously with Wcons', the combination Qmh^ is kept fixed by also varying h. The inclusion of an additional 
tight Gaussian prior on the Hubble parameter {h — 0.72 ± 0.02) breaks the main degeneracy, resulting in the 
tighter inner contours. However, there remains a secondary degeneracy between the value of uicons and the physical 
baryon density Qhh^ ■ This latter quantity can be varied to change the size of the sound horizon and compensate 
for changes in Da{z — 1100) due to Wcons, keeping the angular peak locations fixed. The result is that, even if 
Qm and h axe known precisely, Planck data cannot determine Wcons to a precision better than ±0.1. 



6 



energy at the required accuracy (see Section 
0). 

• Redshift surveys of 10^ — 10^ galaxies in 
the local Universe provide an independent 
measure of via the shape of the galaxy 
power spectrum (if Jigcaiar, ^^b, the neutrino 
masses and any running spectral index are 
known). We argue that the efficacy of this 
cosmological probe is fundamentally limited 
by uncertainties regarding galaxy bias. In 
any case, Figure ^ demonstrates that the 
CMB cannot produce a very accurate mea- 
surement of dark energy even with the ad- 
dition of extra data fixing and/or h. 

In our opinion, the critical additional informa- 
tion required to answer the profound questions 
posed in this Section will be provided by experi- 
ments measuring: 

• Acoustic oscillatory features in the cluster- 
ing power spectrum (Experiment I, Section 

• Weak gravitational lensing by large-scale 
structure, otherwise known as 'cosmic 
shear' (Experiment II, Section 

In this paper we argue that radio surveys with 
the SKA will provide the ultimate databases for 
these experiments, assuming that the telescope is 
designed with a wide enough field-of-view. 

2. Assumptions of this paper 

2.1. SKA capabilities 

We adopt in this paper the SKA capabili- 
ties given by the current 'Strawman Science Re- 
quirements' (Jones 2004) with one key exception: 
the instantaneous (1.4 GHz) field-of-view (FOV) 
must exceed the Strawman value {FOV ^ 1 deg^) 
by more than an order of magnitude for our pro- 
posed experiments. 

AbdaUa & Rawlings (2004; hereafter AR04) 
have considered this requirement in more detail, 
concluding that FOV and f3 (the ratio of the in- 
stantaneous bandwidth to the frequency range 
0.5 — 1.4 GHz required by Experiment I) must 
obey the relation (/3/G.25) x FOV > 40 deg^ for 



the survey to have a reasonable ('^ 1 yr) duration. 
We emphasise that FOV is the field-of-view ob- 
tained at 1.4 GHz (HI at z = 0) and that the cal- 
culations of AR04 assume that this field-of-view 
scales with redshift (oc i^^^ for the purposes of 
this study, where v is the frequency of redshifted 
HI) . AR04 also note the necessity for at least half 
the collecting area to reside in a core of diame- 
ter ^ 5 km for Experiment I, and that Experi- 
ment II requires significant collecting area in the 
longer baselines to deliver the sharp angular res- 
olution (at 1.4 GHz) necessary for measuring cos- 
mic shear. Both of these requirements are already 
in the 'Strawman' design. 

Using these SKA specifications, and assuming 
that the field of view (in units of deg^) scales as 
jy^^, the effective sensitivity of a 'tiled' survey^ is 
sufficient to detect an M^j « 6 x 10^ Mq galaxy 
with a signal-to-noise ratio of 5 out to redshift 
z 1.5 (assuming no evolution in the break of 
the HI mass function; see Sec. 2 of Rawlings et 
al., this volume). 

In this paper we have not attempted to model 
those aspects of our proposed surveys which will 
be very dependent on the precise SKA design, 
such as how the usable FOV scales with wave- 
length for Experiment I, or what dynamic range 
is achievable for Experiment II; we note that in 
practice these limitations might be extremely sig- 
nificant. We will also ignore restrictions on 'all- 
sky' surveys due to finite processing resources 
which - even given likely technological advances 
over the coming decade - may still prove impor- 
tant, given the huge datasets that will be gener- 
ated. 

2.2. Other assumptions 

For the remainder of this paper we assume, un- 
less stated otherwise, a fiducial model consisting 
of a spatially-flat {flk — 0) Universe with cosmo- 
logical parameters h — 0.7, Jim — 0.3, Q^/n^ = 
0.15, (78 = 1 and nscaiar — 1 (broadly consis- 
tent with the recent WMAP results, Spergel et 
al. 2003). We adopt a fiducial dark energy model 
characterized by an equation-of-state w{z) = 

^In which pointings just overlap at 1.4 GHz, and each 
receives a 4-hour integration so that the effective exposure 
time scales as {v/lA GHz)~^. 
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Wq + WiZ where {wq.wi) = (—0.9,0), in order 
to illustrate the accuracy with which the SKA 
will be able to discriminate a general dark energy 
model from a cosmological constant model, for 
which {wq.wi) = (—1,0). Unless stated other- 
wise, we neglect contributions from tensors (pri- 
mordial gravitational waves) and massive neutri- 
nos. We assume that structure in the Universe is 
seeded by adiabatic Gaussian fluctuations."^ 

3. Experiment I: HI emission line surveys 

The SKA will enable revolutionary progress in 
the field of large-scale structure. In particular, it 
will map out the cosmic distribution of neutral 
hydrogen by detecting the HI 21cm transition at 
cosmological distances that are almost entirely in- 
accessible to current instrumentation. Once an 
HI emission galaxy has been located on the sky, 
the observed wavelength of the emission line au- 
tomatically provides an accurate redshift, locat- 
ing the object's position in the three-dimensional 
cosmic web. 

We simulated the co-moving space densities of 
Hl-emitting galaxies as a function of redshift by 
extrapolating the observed local HI mass function 
using model 'C of AR04. We note that, although 
an integral property of the HI mass function is 
constrained by observations of damped-Lya sys- 
tems, considerable uncertainty remains concern- 
ing the redshift-evolution of its break. However, 
reasonable models indicate that the SKA, if de- 
signed with a sufficiently wide field-of-view, can 
survey the entire visible sky in a year of opera- 
tion, locating ~ 10^ HI emission galaxies over a 
volume stretching to redshift z « 1.5 with suf- 
ficient number density that clustering statistics 
are limited by cosmic variance, not by shot noise. 
In this case, the precision with which the galaxy 
clustering pattern can be measured is limited en- 
tirely by the cosmic volume mapped: the SKA be- 
comes the premier large-scale structure machine. 
In comparison, planned next-generation optical 
redshift surveys (e.g. the proposed KAOS spec- 

'^An SKA redshift survey would provide a useful probe of 
non-Gaussianity by detecting large numbers of superclus- 
ters, which correspond to high-peak fluctuations in the pri- 
mordial density field evolving into the quasi-linear regime. 



trograph; http : / /www . noao . edu/kaos) will only 
cover ^ 1000 deg"^ at 2; ^ 1, whereas all-sky 
photometric redshift surveys (e.g. Pan-STARRS; 
http : //panstarrs . if a.hawaii . edu) inevitably 
result in a serious degradation of radial power 
spectrum modes due to the relatively large red- 
shift error.* 

3.1. Galaxy power spectrum 

The resulting HI redshift survey yields an im- 
provement in power spectrum precision of be- 
tween one and two orders of magnitude compared 
to the current state-of-the-art (Figure delin- 
eating features such as the acoustic oscillations 
(see Section In addition, the SKA clustering 
map traces clustering modes on large scales be- 
yond the power spectrum turnover (i.e. on scales 
k < 0.02/1 Mpc~^) that are entirely inaccessi- 
ble to local redshift surveys; together with extra 
small-scale modes (to k w 0.2 h Mpc^^) that are 
hidden by non-linear clustering at redshift zero. 

What cosmological information can we glean? 
The galaxy clustering pattern encodes a wealth 
of data about the underlying constituents of the 
Universe and the basic gravitational processes 
governing the formation of structure. In partic- 
ular, the shape of the power spectrum depends 
on a combination of cosmological parameters. In 
Figure |31 we demonstrate that parameter con- 
straints resulting from the SKA power spectrum 
vastly improve upon those obtained using the 
Sloan Digital Sky Survey (SDSS). Note in partic- 
ular that, since the acoustic oscillations are de- 
tected in the SKA power spectrum (see Section 
I3.2|l . the baryon fraction is strongly constrained, 
whereas the SDSS contours are compatible with 
zero baryons. 

The precise determination of the power spec- 
trum on large scales can be compared with the 
predictions of inflationary models, testing for ef- 
fects such as departures of the primordial power 
spectrum Pprim(fc) from a pure power-law. The 
narrow window function in fc-space (resulting 

*Even with a very optimistic redshift error Az = 0.02, 
a photometric redshift survey must cover ^ 10 times the 
area of a spectroscopic experiment to map out the same 
number of Fourier modes in the linear regime (Blake &; 
Glazobrook 2003). 
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Figure 2. Simulated measurement of the galaxy clustering power spectrum by a 1-year SKA survey, compared 
to the 2dFGRS 100k data release (Percival et al. 2001). The measurements have been normalized to the same 
underlying power spectrum model, keeping the (la) fractional errors in each fc-bin the same. For the SKA data, 
we assume that the redshift evolution of clustering may be divided out so that results from different redshift slices 
can be combined. We note that the 2dFGRS P{k) data points are heavily correlated in this binning, whereas 
the SKA survey possesses a much narrower 'fc-space window function' owing to the vastly greater cosmic volume 
probed, implying that the resulting P{k) measurements are barely correlated. The SKA clustering data delineates 
the power spectrum turnover on large scales {k < 0.02 h Mpc~^) and the acoustic oscillations (Section 13.21 . as a 
direct result of mapping a vastly greater cosmic volume V. The error in a power spectrum measurement in the 
cosmic- variance- limited regime scales as l/\/V, and Vska ~ 500 V2dF. 
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Figure 3. Joint constraints on the cosmological parameters (f2mft, /b ~ i^h/^m) resulting from fitting linear 
theory models (calculated using the CAMB software; Lewis, Challinor & Lasenby 2000) to the shape of the SKA 
galaxy power spectrum (assuming a linear galaxy biasing scheme). Contours are shown for the analyzed Sloan 
Digital Sky Survey (SDSS) data (Pope et al. 2004; la, 2a, 3a) and for a simulated SKA HI emission-line survey 
(inner contour; 2a). We only use SKA clustering data in the linear regime (with k < 0.2 h Mpc~^). In order to 
match the assumptions of the SDSS contours, we fix riscaiar ~ 1 and marginalize over only the Hubble constant. 
Whilst the SKA produces an exquisite delineation of the power spectrum, we argue that this Figure may be wildly 
optimistic because the bias model cannot be assumed to be linear at the required level of precision. 
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from the vast cosmic volume probed by the SKA) 
imphes that any sharp features in Pprim(fc) Stve 
not smoothed out and may be detected in a man- 
ner not possible using the CMB. It has been sug- 
gested that the intriguing outliers in the CMB 
temperature power spectrum measured by the 
WMAP satellite (at ^ - 30 and € ~ 200) could 
be due to unknown physics on scales smaller than 
the Planck length, which are imprinted during 
inflation (e.g. Martin & Ringeval 2004). The 
SKA would delineate any such effects. Further- 
more, comparison with large-scale CMB modes 
can yield information about any tensor contribu- 
tion to the CMB anisotropics, arising from grav- 
ity waves in the early Universe.^ This is a sec- 
ond potential means, together with pulsar obser- 
vations, by which the SKA can be used as a grav- 
itational wave detector. 

However, we wish to argue that the appli- 
cations of the SKA large-scale structure mea- 
surements discussed above are potentially flawed 
(and at best severely complicated) by our lack of 
knowledge of the detailed biasing scheme linking 
the observed galaxy distribution to the underly- 
ing pattern of dark matter fluctuations produced 
by theoretical models (e.g. Coles 2004). In gen- 
eral, galaxy biasing is almost certainly stochastic, 
scale-dependent, non-local and non- linear (Dekel 
& Lahav 1999): the assumption of a simple lin- 
ear bias, implicit in Figure|31 is likely to be inade- 
quate when confronted with sub-percent precision 
power spectrum data. 

3.2. Acoustic oscillations 

Given the likely complexities of galaxy bias, we 
seek cosmological experiments that are less sensi- 
tive to the biasing model. The power spectrum of 
galaxies on large scales ( > 30 Mpc) should con- 
tain a series of small-amplitude acoustic oscilla- 
tions of identical physical origin to those seen 
in the CMB.^° These features result from oscil- 



^Assuming that such a tensor contribution may be disen- 
tangled from any large-scale effects due to clustering of 
dark energy, galaxy bias, and the potential presence of 
'isocurvature' density fluctuations in the early Universe 
(in addition to the usual adiabatic perturbations). 
-'^''The acoustic oscillations imprinted in the matter power 
spectrum possess significantly lower amplitude than those 
found in the CMB power spectrum, because the baryonic 



lations of the photon-baryon fluid before recom- 
bination, and encode a characteristic scale - the 
sound horizon at recombination - accurately de- 
termined by CMB observations as described in 
Section^ This scale can act as a standard ruler 
(Eisenstein 2002). Its recovery from a galaxy red- 
shift survey depends on the assumed cosmological 
parameters, particularly the dark energy model, 
and thus constrains that model over a range of 
redshifts (Blake & Glazebrook 2003; Seo & Eisen- 
stein 2003). 

The application of this cosmological test does 
not depend on the overall shape of the power 
spectrum, which can be divided out using a 
smooth fit, but only on the residual oscilla- 
tory signature of the acoustic peaks. Hence the 
method is insensitive to smooth broad-band tilts 
in P{k) induced by such effects as a running spec- 
tral index, redshift-space distortions, and com- 
plex biasing schemes (e.g. based on haloes; Seljak 
2000). It would be very surprising if the biasing 
scheme introduced oscillatory features in /c-space 
liable to obscure the distinctive acoustic peaks 
and troughs. In the absence of major system- 
atic effects, constraints on the dark energy model 
are limited almost entirely by how much cosmic 
volume one can survey, rendering this an ideal 
experiment for the SKA, which can map the HI 
galaxy distribution out to z « 1.5 (where three 
acoustic peaks lie in the linear regime). 

In Figure^we display simulated power spectra 
for a survey covering 20,000 deg^, divided into 
redshift slices. If the observer applies the wrong 
dark energy model when converting HI redshifts 
to cosmic distances, then the peaks and troughs 
in the recovered P{k) lie in the wrong places: the 
experiment becomes a powerful probe of dark en- 
ergy. 

The standard ruler provided by the acoustic 
oscillations may be applied separately to the tan- 
gential and radial components of the power spec- 
trum, resulting in independent measurements of 
the co-ordinate distance x{z) to each redshift slice 
and its rate of change x'{z) = dx/dz, respec- 



component of the matter distribution, which contains the 
oscillatory imprint, is sub-dominant to the CDM compo- 
nent. 
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Figure 4. Simulated galaxy power spectra for an SKA survey of 20,000 deg^, analyzed in redshift slices of width 
Az = 0.2. Each power spectrum is divided by a smooth polynomial fit, revealing the sinusoidal imprint of acoustic 
oscillations, and is shifted along the j/-axis for clarity. As redshift increases, the linear regime extends to smaller 
physical scales (larger values of k), unveiling more peaks. The acoustic oscillation 'wavelength' may be used as 
a standard ruler. This is illustrated by the Figure: the true dark energy model is w — —1, but the observer 
has incorrectly assumed w = —0.8 when constructing the power spectra, and consequently the recovered acoustic 
scale disagrees with that observed in the CMB (represented by the vertical dotted lines). 
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tively.^^ We fitted these simulated measurements 
of x{z) and x'{z) over a series of redshift slices, 
assuming a dark energy model w{z) = wq + wiz. 
Figure \^ illustrates the accuracy of recovery of 
(wq, wi) assuming a fiducial cosmology (—0.9, 0), 
demonstrating that in this case we can dismiss 
a cosmological constant model (—1,0) with high 
significance. The accuracy of our prior knowl- 
edge of il„i and h affects the tightness of the dark 
energy constraints. These parameters must be 
known to standard deviations a{flm) = 0.01 and 
a{h) = 0.01 if these uncertainties arc not to domi- 
nate the overall error in {wq, wi). These demands 
are not unrealistic on the SKA timescale: SKA 
observations of masers in the vicinity of black 
holes will determine h to ^ 1% accuracy (Green- 
hill, these proceedings) and CMB data from the 
Planck satellite will fix Vlmh^ to a similar preci- 
sion (Balbi et al. 2003). 

Given that the SKA will be unavailable until 
~ 2020, what other projects plan to survey the 
acoustic peaks? Two leading contenders are the 
multi-object near- infrared/optical spectrographs 
FMOS and KAOS, due to commence operation 
later this decade at the Subaru and Gemini tele- 
scopes, respectively.^^ These projects may per- 
form redshift surveys of a few hundred deg^, de- 
riving constraints on {wq,wi) that are several 
times poorer than the 20, 000 deg^ SKA survey 
presented above (the error bars roughly scale as 
l/^Z/sky)- In radio wavebands, SKA prototypes 
may be developed such as the HYFAR proposal 
(Bunton et al. 2003). The speed of surveying 
= 1 HI emission galaxies for such a prototype 
would lag the SKA by roughly an order of magni- 
tude; furthermore, the bandwidth may be signif- 
icantly smaller than that envisaged for the SKA. 

^^The actual quantities determined by the experiment are 
the ratio of x and x' to the sound horizon s at recombina- 
tion (which determines the characteristic spatial scale of 
the acoustic oscillations). The value of s is determined by 
CMB observations, via Qrah? and Q,\^h? , but is affected to 
some extent by uncertainities in Cm and h (Figure|SJ. The 
tangential component of P{k) measured in a survey slice 
at redshift z determines x{z) because this latter quantity 
governs the tangential scale of the recovered galaxy distri- 
bution in that redshift slice: displacements Ar are deter- 
mined by Ar = xl\d. The value of x'{z) fixes the radial 
scale of the recovered distribution; Ar = (dx/dz) Az. 
i^KAOS is still unfunded. 



In summary: even if these precursor experi- 
ments succeed in detecting a deviation from Ein- 
stein's cosmological constant, the SKA will still 
be required for a precise delineation of the new 
model. The key advantage offered by a radio tele- 
scope is that it may be designed with a field of 
view exceeding an optical spectrograph by a fac- 
tor of 100: the SKA can survey all the available 
volume by studying the whole visible sky. More- 
over, SKA sensitivity is necessary and sufficient 
for detecting tracers of this volume out to z = 1.5. 

4. Experiment II: radio continuum surveys 
for cosmic shear 

Light rays follow geodesies, which bend in the 
presence of matter. It follows that a coherent 
shape distortion is imprinted in the distribution 
of distant background galaxies by mass fluctu- 
ations in the intervening cosmic web: there is 
a small tendency for galaxies to be tangentially 
aligned. This pattern of cosmic shear encodes a 
vast body of information and, in principle, is a re- 
markably powerful cosmological probe. The shear 
pattern is directly sensitive to the mass distribu- 
tion predicted by the underlying theory, and docs 
not depend on the complex details of how galaxy 
light traces mass. Furthermore, cosmic shear may 
be used to measure dark energy: which controls 
both the angular diameter distances intrinsic to 
the lens equation (e.g. the source-lens distance) 
and the redshift evolution of the mass fluctuations 
producing the lensing. 

A successful cosmic shear survey has three lead- 
ing design requirements: very high image quality 
for reliable shape measurements, a high source 
surface density ( > 100 arcmin"^) to limit sta- 
tistical noise, and a wide survey area to reduce 
cosmic variance. Each of these demands is deliv- 
ered superbly by an SKA radio continuum survey 
(Schneider 1999). 

The leading systematic for ground-based op- 
tical cosmic shear experiments is the difficulty 
in quantifying systematic variations in the point- 
spread function, due to inevitable changes in the 
atmospheric seeing and telescope properties with 
position and time. The psf determination is ac- 
complished by observing stars in the field (of 
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Figure 5. Constraints on a dark energy model w(z) — wo + wiz achievable by an SKA survey of the acoustic 
peaks, illustrated by sets of (68%, 95%) contours corresponding to different assumed priors. The solid contours 
correspond to perfect knowledge flm = 0.3 and h — 0.7. For the dashed contours we impose Gaussian priors 
with standard deviations a{^lm) = 0.01 and (t(/i) — 0.01, which are not unrealistic in the SKA epoch. In each 
of these two cases we are able to distinguish a fiducial model {wq,Wi) — (—0.9,0) from a cosmological constant 
(marked by the solid circle). Weaker Gaussian priors a(i^m) = 0.03 and o-(/i) = 0.05 (the dot-dashed contours) 
are inadequate for this purpose. Note that there remains a residual degeneracy in the («;o,wi) plane because 
at a given redshift, approximately the same cosmology can be obtained by increasing wq and decreasing wi. 
This degeneracy is partially broken by measurements across a range of redshift slices (and by the independent 
determinations of x and x' for each redshift). 
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which there are a hmited number density). In 
comparison, the point-spread function of a ra- 
dio telescope is well-determined and stable (being 
simply derived from the interferometer baseline 
distribution, i.e. the synthesized beam). In addi- 
tion, the planned SKA angular resolution (« 0.05 
arcsec at 1.4 GHz) vastly improves upon that ob- 
tainable over wide fields from the ground. 

Space-based optical cosmic shear surveys 
should also exist by the time the SKA is avail- 
able. However, the planned area for these surveys 
is « 300 deg^ (a SNAP-like experiment; Massey 
et al. 2004), which will be dwarfed by an SKA all- 
sky continuum survey. The fraction of sky /sky 
mapped is directly related to the precision of the 
shear power spectrum measurement via a factor 

1/ \/ /sky 

An SKA continuum survey will achieve sensi- 
tivities of ~ 30 njy in a 4-hour pointing, far ex- 
ceeding depths attained by contemporary radio 
surveys. In order to simulate the results of an 
SKA cosmic shear experiment, we must assume 
a model for the radio source populations at these 
unprecedented flux densities. We generated this 
model from a combination of two populations: 
starburst galaxies and quiescent spiral disks (we 
ignore for now the Active Galactic Nuclei which 
dominate at the highest flux densities but which 
have too low surface density to contribute signif- 
icantly to the weak lensing signal at SKA sen- 
sitivities). We extrapolated local radio luminos- 
ity functions for these two populations (Sadler et 
al. 2002) as a function of redshift using a Press- 
Schechter-based approach, and folded in a pre- 
scription for the physical sizes of each population. 
Full details of these models will be published else- 
where (Abdalla et al., in prep.). 

Using these population models, we simulated 
radio skies as a function of SKA integration time 
and angular resolution. Following the method 
of Massey et al. (2004), we then employed stan- 
dard shape-estimation software (Kaiser, Squires 
& Broadhurst 1995) to recover the shear disper- 
sion per galaxy arising from intrinsic shapes and 
from measurement errors (a-y) together with the 
surface density of 'usable galaxies' (ug) after cuts 
for objects that are unresolved or faint. The val- 
ues achieved, ~ 0.2 and ng ^ 500 arcmin"^. 



are comparable to space-based optical surveys 
(Massey et al. 2004) whilst covering a survey area 
100 times larger. Figure|n|illustrates the measure- 
ment accuracy of the cosmic shear angular power 
spectrum using these data.^^ 

The resulting measurements of the dark energy 
model can be estimated using the Fisher matrix 
methodology (e.g. Refregier et al. 2004). Figured 
illustrates joint constraints on {wo,wi) using the 
shear power spectrum up to ^ = 10^. The anal- 
ysis requires models for both the source redshift 
distribution (which can be readily inferred from a 
deep HI emission-line survey of a sub-area of the 
sky) and the linear theory mass power spectrum 
(e.g. calculated using CMBFAST; Seljak & Zal- 
darriaga 1996) corrected for non-linear effects in 
accordance with the prescription of e.g. Smith et 
al. (2003). 

The power of this experiment may be increased 
by splitting the lensed galaxies into two broad 
redshift bins (Refregier et al. 2004). This is 
readily, if roughly, accomplished by separating 
the radio continuum sources into two groups de- 
termined by detection (or absence) in the HI 
emission-line survey of Experiment I. 

We note that a cosmic shear survey does not 
necessarily require exquisitely-measured galaxy 
shapes, because there is an irreducible scatter 
resulting from the unknown ellipticity of each 
galaxy before shear. Instead, the 0.05-arcsec 
SKA resolution is required to (1) resolve starburst 
galaxy disks at high redshift and (2) eliminate 
radio source confusion. A potential residual sys- 
tematic error may lie in intrinsic alignments of 
galaxy disks (owing to structure formation pro- 
cesses) which may masquerade as coherent shear 
signal. For example, some starburst galaxies may 
exist in pairs, forming merging disk systems. Tak- 
ing an extremely cautious approach, we assume 
for now that starburst galaxies would not be used 
in the cosmic shear analysis. This reduces the sur- 
face density by a factor of two, but the SKA still 
yields a vast improvement in dark energy con- 



^^We note that it is also possible to analyze a radio cosmic 
shear survey directly in the uv-pla.ne (Chang & Refregier 
2002). 

^^In the ideal cosmic shear experiment, the unsheared ob- 
jects would all be circular. 
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Figure 6. The cosmic shear angular power spectrum measured from a simulated SKA radio continuum survey. 
The upper (red) curve is for a model with Sim — 0.25 and Wcons = —0.9; the lower (blue) curve is for = 0.3 
and Wcons = —1; the dotted lines show the ±1(t range for f2m = 0.3 and Wcons = —0.9. We assume that the survey 
contains a surface density of usable galaxies Ug — 200 arcmin"^ and covers one hemisphere (/sky ~ 0.5). 
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Figure 7. (68%, 95%, 99%) likelihood contours for a dark energy model ■w{z) = wo + wiz from an SKA cosmic 
shear survey over a hemisphere, assuming 200 usable galaxies per square arcminute. We marginalize over the 
other cosmological parameters (f2m, fib, h, Us). These contours are conservative because no redshift information 
per galaxy has been included. 
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straints compared to future space-based surveys. 

The shear angular power spectrum is the 
most basic lensing experiment; we note that 
the bi-spectrum (three-point statistic) promises 
comparable, statistically-independent constrain- 
ing power on small scales (Takada & Jain 2004). 
A third independent approach is to target the 
most massive clusters around which the imprint 
of cosmic shear is most significant. Deep HI 
follow-up can yield full background source red- 
shift information, allowing us to cross-correlate 
the cosmic shear signal in redshift slices be- 
hind the massive lens, yielding a purely geomet- 
ric probe of cosmic distances ('cross-correlation 
cosmography', e.g. Jain & Taylor 2003). This 
method provides a probe of the angular diame- 
ter distance as a function of redshift independent 
of the acoustic oscillations in the clustering power 
spectrum. However, a potential systematic error 
is introduced by the presence of additional pro- 
jected masses along the line-of-sight. 

5. The SKA as the Dark Energy Machine 

The current frontier of cosmological research is 
the characterization of the properties of the 'dark 
energy' that dominates the energy density of to- 
day's Universe and drives its accelerating expan- 
sion. In this paper we have focussed on two key 
experiments - acoustic oscillations (Experiment 
I) and cosmic shear (Experiment II) - by which a 
wide-field SKA can generate measurements of the 
dark energy model that are significantly more ac- 
curate than any competition by ~ 2020. Even if 
precursor surveys detect evidence for a deviation 
from Einstein's cosmological constant, the SKA 
will be required for a precise delineation of the 
new model. 

In Figure IHl we compare overall SKA measure- 
ments of a dark energy model w{z) = wq + wiz 
with (arguably) the most promising competitor: 
the (currently unfunded) satellite mission SNAP, 
planned to laimch in ^ 2014, whose basic goal is 
to locate approximately two thousand supernovae 
to z K, 1.7 ( |http : // snap . Ib l . gov| Aldering et 
al. 2004). 15 

^^We note that the planned SNAP mission also includes a 
weak lensing component that will yield outstanding space- 



Observations of distant Type la supernovae as 
'standard candles' provided the first clear evi- 
dence for the accelerating expansion of the Uni- 
verse (Riess et al. 1998; Perlmutter et al. 1999). 
Future high-redshift supernova searches will be a 
powerful probe of the dark energy model. How- 
ever, we argue that it is not yet proven that the 
dimming effect of cosmic acceleration can be dis- 
tinguished from other systematic effects with suf- 
ficient precision to render this cosmological probe 
as reliable as Experiments I and II above. Such 
systematic effects include changes in the intrin- 
sic properties of supernovae with galactic envi- 
ronment (e.g. metallicity), dust extinction from 
both the host galaxy and the Milky Way^^, and 
the difficulties of photometric calibration (super- 
nova /-C-corrections, together with the inadequacy 
of local 'standard stars' such as Vega for provid- 
ing sub-percent-level flux calibration) . 

Many other probes of dark energy have been 
suggested in the literature (see e.g. Gurvits, this 
volume). We argue that, whilst these probes 
will provide important cross-checks of the cosmo- 
logical model, none can yield the raw precision 
or control of systematic errors afforded by our 
Experiments I and II. For example, using SKA 
(or other) data we can perform number counts 
of haloes as a function of circular velocity and 
redshift to infer the cosmological volume element 
(Newman & Davis 2000). But this experiment 
requires an accurate theoretical prediction of the 
number density distribution with circular veloc- 
ity (i.e. a precise understanding of baryons) to- 
gether with precisely-measured rotation curves at 
high redshift. Alternatively, number counts of 
galaxy clusters can be utilized (e.g. Weller, Battye 
& Kneissl 2002). The main systematic difficulty 
here is connecting the cluster observable (e.g. X- 
ray luminosity or SZ decrement) to the halo mass 
(upon which the cluster abundance is exponen- 

based image fidelity. However, as noted in Section |3 an 
SKA cosmic shear survey can cover a area ~ 100 times 
greater with comparable resolution. 

^^The recent measurement of several supernovae at z > 1 
by the Hubble Space Telescope (Riess et al. 2004) has as- 
suaged this doubt to some extent by detecting a transition 
with increasing redshift from accelerating to decelerating 
expansion, which would not occur if the inferred accelera- 
tion was an artefact of dust extinction. 
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Figure 8. Comparison of constraints on the dark energy equation-of-state, parameterized as w{z) — wo + wiz, 
obtained by the proposed satellite mission SNAP (using both the supernova and cosmic shear experiments) and the 
SKA (via acoustic oscillations and cosmic shear). The contours represent 'ilcr' confidence intervals. A Gaussian 
prior upon the matter density, a{Qni) = 0.01, was consistently applied to the supernova and acoustic peaks 
analyses. For the latter we also assume a prior (t(/i) — 0.01 for the Hubble parameter, which could be obtained 
by measuring masers using the SKA. We assumed that the SNAP supernova component produces measurements 
of the luminosity distance to an accuracy of 1.4% in each of 50 redshift bins between z — 0.1 and z = 2. The 
cosmic shear contours for both SNAP and SKA are conservative because the sample has not been split into two 
or more bands in redshift. This analysis method was shown by Refregier et al. (2003) to improve significantly the 
constraints from SNAP cosmic shear, and would have a similar effect on the SKA measurements. 
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tially sensitive). The frequency of strong gravita- 
tional lensing (the number of lensing arcs in the 
vicinity of massive haloes) is also sensitive to dark 
energy (Bartelmann et al. 2003), but in manner 
very dependent on accurate simulations (see also 
Koopmans, Browne & Jackson, this volume)DW. 
We argue that these various experiments will pro- 
vide insight into the complex processes of struc- 
ture formation, rather than into the underlying 
cosmological model. 

Finally, if matter does not dominate the dy- 
namics of the Universe, then gravitational poten- 
tial wells begin to 'decay', inducing a net shift 
in the energies of CMB photons passing through 
these wells. This phenomenon is known as the 
late-time Integrated Sachs Wolfe (ISW) effect. 
Its influence on the CMB power spectrum is re- 
stricted to large scales {£ < 10) where it is clouded 
by cosmic variance. However, its presence induces 
a correlation between CMB anisotropics and the 
low-redshift matter distribution, which has been 
detected using current radio surveys (e.g. Boughn 
& Crittenden 2002) and could be probed further 
with the SKA. Measurement of the late ISW ef- 
fect is an important confirmation of the dark en- 
ergy model, and is sensitive not just to the evolu- 
tion of dark energy, but also to its clustering prop- 
erties (e.g. Weller & Lewis 2003; Bean & Dore 
2003). 

We conclude that, owing to its capacity to 
undertake cosmic surveys on an unprecedented 
scale, the SKA will become the paramount tool 
for answering the key questions of the era of preci- 
sion cosmology. The SKA can perform clean and 
powerful experiments to characterize dark energy, 
mapping the acoustic oscillations in the cluster- 
ing power spectrum and the patterns of cosmic 
shear distortion in the radio sky. Furthermore, 
SKA data can enlighten other pressing cosmolog- 
ical questions, pursuing tests of inflationary theo- 
ries and measuring spatial flatness with unprece- 
dented precision. 
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